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k*- ' Abstract 
CO 

A modified Friedmann equation which arises in an extension of general relativity 
which accommodates a time-dependent fundamental length L(t) leads to cosmolog- 
ical models where the scale factor diverges with an essential singularity at a finite 
future time. Such models have no dark energy in the conventional sense of energy 
possessing a truly simple pressure-energy relationship. Data on supernovae restrict 
the time from the present until the Rip to be generically longer than the current age 
of the Universe. 
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1 Introduction 



In this talk I report on the recent joint work with Takahashi[l . 

The cosmic concordance of data from three disparate sources: Cosmic Microwave Back- 
ground (CMB), Large Scale Structure (LSS) and Type la Supernovae (SNela) suggests 
that the present values of the dark energy and matter components, in terms of the critical 
density, are approximately fix — 0.7 and fi m ~ 0.3. 

The equation of state of the dark energy w = p/p suggests the possibility that w < — 1, 
first studied by Caldwell|2j and subsequently in a number of papersjSl Ej. 

The conclusion about the make-up of our Universe depends on assuming that general 
relativity(GR) is applicable at the largest cosmological scales. Although there is good 
evidence for GR at Solar-System scales [5 there is no independent evidence for GR at 
scales comparable to the radius of the visible Universe. The expansion rate of the Universe, 
including the present accelerating rate of cosmic expansion can be parameterized in the 
right-hand side of the Friedmann equation by including a dark energy density term with 
some assumed time dependence on the scale parameter a(t): pde ~ oP with j3 = — 3(l+w). 
This is a rather restricted function if we assume that the equation of state w is time- 
independent. But as soon as we admit that it may depend on time w(t) then the function 
on the right hand side of the Friedmann equation becomes completely arbitrary, just as 
does its designation dark energy density. 

Present data are fully consistent with constant w = —1 corresponding to a cosmological 
constant which we may accommodate on the left hand side of the Friedmann equation 
describing the expansion rate of the Universe. But cases with to^-l, including w < — 1, 
are still permitted by observations. In this case, there is a choice between cooking up a 
"dark energy" density with a particular time dependence pde ~ on the right-hand 
side of the Friedmann equation or changing the left-hand-side by changing the relationship 
between the geometry and the matter density, i.e. by changing GR. Even if we do the 
latter, from the viewpoint of the Friedmann equation, we can always find a time-dependent 
term on the right-hand-side which is equivalent and which we may call "dark energy" 
and thus preserve GR in some form. Eventually, this distinction may come down to 
observational tests of whether a particular change in the geometry predicted by GR can 
be detected, other than by the expansion rate of the Universe. 

The case constant w < — 1 has the interesting outcome for the future of the Universe 
that it will end in a finite time at a "Big Rip" before which all structure disintegrates [6 . 

In the present article, we shall study an amalgam of the modification of GR due to 
Dvali, Gabadadze and Porrati[7j (DGP) and the idea of a Big Rip, in fact here a Bigger 
Rip. This will be based on admittedly ad hoc ansatz for terms in modified Friedmann 
equations but the results are sufficiently interesting to examine and such modifications 
may be constrained by observational data. 
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2 Set up 



The DGP gravityjZ] arises from considering the four- dimensional gravity which arises from 
five-dimensional general relativity confined to a brane with three space dimensions. The 
underlying action is: 

S = M{tf J d 5 xVGK® + M 2 Planck J d'x^gR (1) 

where and R are the scalar curvature in 5- and 4- dimensional spacetime respectively, 
and G and g are the determinant of the 5- and 4- dimensional spacetime metric. 

This leads to an interesting modification of GR which embodies a time-dependent 
length scale L(t) = Mp lanck /M(t) 3 . For cosmology it is natural to identify, within a 
coefficient of order one, L(t) with the Hubble length L[t) = H(i)~ l at any cosmological 
time t. Actually we are slightly generalizing the original DGP approach to include time 
dependence of the length scale L(t). 

Taking the four dimensional coordinates to be labeled by i, k = 0, 1, 2, 3 leads to the 
following modification of Einstein's equation: 



#k_l Rg ik\ + V(/ 



2 J L(t)y/g LV 2 



(2) 



where [(...)] means: 



dx[(f(x))] = f'(0)S(x). (3) 



It is interesting to generalize the Schwarzschild solution to this modification of GRjHJ. 
One finds that the modification of the Newton potential at short distances is given by: 



V{r) 



Gm 4 V Gmy/r 



2r 



where r g = 2Gm is the Schwarzschild radius. 



L(t) 
L(t) 



(4) 



The fractional change in the Newtonian gravitational potential at cosmological time t 
at orbital distance r from an object with Schwarzschild radius r g is therefore 



AV 



V 



L(t)\ 



(5) 



In the Bigger Rip scenario we will describe the characteristic length L(t) will decrease 
with time according to 

nt) = Kt^ntf (6) 
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where the power satisfies p > 1 (p < 1 implies that L(t) would increase) and where 

T(t) = ( 7 ) 

in which t r i p is the time of the Rip. A bound system will become unbound at a time tu 
when the correction to the Newtonian potential becomes large. We make adopt the value 
of tu defined from Eq.© by 

8r 3 



L{tu) 2 r 9 



(t rip -tu) = - (9) 



We can rewrite Eq.(jSJ) as: 

(t r i p - tu) = x ., 

where 7 = (t rip - to) -1 - 

We shall define another later time t caus as the time after which the two objects of a 
bound system become causally disconnected from t caus until t r % v . This is defined by the 
equation: 

(t r ip t caus ) I — r J (10) 



c \ a(ti 



As an example taking p = 1 with the values Lq = H 1 = (lAGy) 1 = 1.3 x 10 28 cm 
and 7 = (20Gy)^ 1 we arrive at the entries in the following Table: 



Bound system 


l {cm) 


r g (cm) 


(trip tu) 


(t r ip t caus ) 


Typical galaxy 


5 x 10 22 


3 x 10 16 


lOOMy 


4My 


Sun-Earth 


1.5 x 10 13 


2.95 x 10 5 


2mos 


31hr 


Earth-Moon 


3.5 x 10 1U 


0.866 


2weeks 


lhr 



Note that the values we find for (t rip — tu) are consistent with those found in [6]. The 
corresponding dark energy would have equation of state w = —1 — I7A) = —1-466 which, 
like that of [B], is now outside of the range allowed by observations [9 J if we assume a 
constant equation of state although it is allowed in the present model with its time- 
dependence. As another example, with more normal present w, we can increase the time 
to the Rip to 7 = (50G?/) -1 in which case w(t ) = —1.19 and the Table is modified to: 



Bound system 


l (cm) 


r 9 (cm) 


(trip tu) 


(trip t caus ) 


Typical galaxy 


5 x 10 22 


3 x 10 16 


250My 


7My 


Sun-Earth 


1.5 x 10 13 


2.95 x 10 5 


5mos 


2days 


Earth-Moon 


3.5 x 10 10 


0.866 


lmo 


2hrs 
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so with the more lengthy wait until the Big Rip the disintegration of structure and causal 
disconnection occur correspondingly earlier before the eventual Rip. 



3 The Bigger Rip 



The modified Friedmann equation for DGP gravity is 
so that we arrive at: 

i = H(t) = H(to)± (12) 

In Eqs. (jllll2|) we can neglect, for the future evolution, the term (p M + p 7 ) / (3Mp lanck ) 
on the right-hand-side of the modified Friedmann equation. Defining 7 = —dT/dt = 
(trip - h)' 1 gives: 

,T(t) dT 



^=-7, ^mfi < 13 > 

and hence, for p — 1, which is similar to dark energy with a constant w < — 1 equation of 
state: 

a(t) = T~^k) (14) 
while for the Bigger Rip case p > 1 one finds 



a(t) = a (t ) exp 



I ^ 1 



rp-i ; (p-i) 7 L(t )J 



(15) 



Here we see that the scale factor diverges more singularly in T for p > 1, hence the 
designation of Bigger Rip. In particular we study the values p — 2, 3, • • • as alternative to 
the "dark energy" case p — 1. 
Inverting Eq. ([T5|) gives: 



T= [l + (p- l)7^o) In a(t)]~T?=TT (16) 

In this case there is strictly no dark energy, certainly not with a constant equation of state, 
but we can mimic it with a fictitious energy density pl by noticing that if 2 ~ T~ 2p and 
writing 

Pl~ [l + (p- l)7^o) In a(t)](^> (17) 
If we use Eqs. (115)1 and (|T7|) in conservation of energy 

1(^3) = _ p | (a3) = ^(^1^ (18) 
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we find a time-dependent wz,(i) for the "fictitious" dark energy 

2 dL(t) , . 

= _1_2 (2Q) 

31 + (p-l) 7 L(t )lna(t) V ; 

so the effective u>z(t) has the limiting values w^it^) = — 1 — |p(7L(i )) and WL(t rip ) = — 1. 

We may check consistency with the space-space components of Einstein's equations 
which are 

§ = -££-4** (21) 

which leads to 

if = -47rG(p L +p L ) = -4irp L (l + w L ) (22) 
so that with H = L^ 1 and Pl = 3H/ (8nGL) we find a wz,(t) consistent with Eq.(j2Hl)- 



Keeping the value L = if _1 = [\AGyY 1 = 1.3 x 10 28 cm and putting 7 = (20Gy) _1 
and p = 2 we arrive at the entries in the following Table: 



Bound system 


l (cm) 


r 9 (cm) 


(trip tjj ) 


(t C aus tjj J 


Typical galaxy 


5 x 10 22 


3 x 10 16 


2.37Gy 


1.14Gy 


Sun-Earth 


1.5 x 10 13 


2.95 x 10 5 


9.6 x W 4 y 


7y- 


Earth-Moon 


3.5 x 10 10 


0.866 


2.5 x 10 4 y 


Qmos. 



Note that for the p = 2 case we have tabulated the difference (t caus — tjj) rather than 
(t r ip — tcaus) because in this case the expansion is so rapid. 



Next we turn to observational constraints on the parameters L and 7 for p = 2. 
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4 Observational Constraints 



In the previous section, we discussed the future universe in the model. In this section, we 
discuss constraints on model parameters from SNelA data. To discuss the constraint, we 
have to include other component such as cold dark matter (CDM) and baryon. Including 
all components, we can write the Friedmann equation as [10 j 



Pm 


1 


3^Planck 


f 4L 2 


= Pm/ Pcrit 


— PmO 




If we define the density parameter fl m = p m j 'p crit = p m o(l + z) 3 , we can rewrite Eq. 



as 

where Qk and Ql are defined as 



H 2 = Hi 



(24) 



■k 



Thus at the present time, we have the relation among the density parameters, 

tt k + (yvT L + y/n L + Qm) 2 = 1. (26) 

To obtain a constraint from the SNela observations, we assume that red-shift depen- 
dence of Ql is as in Eq. (JT7jl and that L is dependent on time as Eq. (|B|); thus we can say 
that we consider a new component p^ which is defined as Eq. (|17|). 

Now we discuss the constraint on this model from SNela data using recent result jH]. 
In Fig. [TJ we show contours of 95 and 99 % C.L. in flL-fl m plane for jLq = (which is 
the constant L case), 0.5 and 1 with p = 2. In the figure, we also plot the line for the 
flat universe. Notice that the line is different from the standard case because we have the 
modified Friedmann equation in this model. To obtain the constraint, we marginalize the 
Hubble parameter dependence by minimizing x 2 f° r the fit. 



In Fig. |2l we show the constraint on Q m - , yLQ plane assuming the flat universe (Qk = 0). 
If we take the value fi m = 0.3, we can find an upper limit for ^L which is 7-L ^ 0.7. This 
implies that the time remaining from now until the Rip is constrained to be generically 
at least somewhat longer than the current age. If we make L larger than the length 
corresponding to the age of the Universe then the upper bound on 7 diminishes and hence 
the time until the Rip increases. 

We note that because the effective equation of state w(t) is varying with time its present 
value w(to) can be more negative than allowed by constraints derived from assuming 
constant w. Our constraint on < 0.7 permits w(t ) = —1 — |p7-L to be as negative 
as w(to) = —1.9 for p — 2. Assuming constant w, on the other hand, gives |«;> —1.2. 
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5 Discussion 



The present article is a natural sequel to our previous paper jl] about dark energy in which 
it was pointed out that no amount of observational data can, by itself, tell us the fate of 
dark energy if we allow for an arbitrarily varying equation of state. The three possibilities 
listed were: there may firstly be a Big Rip, or secondly dark energy may dominate but with 
an infinite lifetime or thirdly the dark energy may eventually disappear leaving a matter- 
dominated Universe. Given that observational data are insufficient, only a successful and 
convincing theory of the past may inform of the future of the Universe. 

The Big Rip was the most exotic of the fates and there seemed tied to a phantom 
w < — 1 dark energy. However, here we have studied a Bigger Rip, in which the scale 
factor is even more divergent at a future finite time than for the Big Rip, which is achieved 
by modifying gravity and omitting dark energy. In the model, as with the phantom 
case, structures become unbound and subsequently their components become causally 
disconnected before the Universe is torn apart in the Rip. 

If we allow an arbitrarily- varying equation of state any new term on the left-hand-side of 
the Friedmann equation can apparently be taken to the right-hand-side and reinterpreted 
as a dark energy. But this will, in general, lead to a relationship between pressure and 
energy which is not truly simple, but very complicated as in Eq. (J2DJ) above. We use the 
term "dark energy" when the time dependence of such a term governing the expansion 
rate of the Universe is sufficiently simple, e.g. it could be required that w(t) is either 
constant or at most polynomial in t. This is not the case for the present model and hence 
justifies the title chosen for this article. 
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Figure 1: Constraint from SNela observation in Q L -Q m plane for 7L0 = (bottom), 
7L0 = 0.5 (middle) and 7L0 = 2 (top). Contours are for 95 % (dotted line) and 99 % 
(dashed line) C.L. constraints respectively. The solid line indicates parameters which give 
a flat universe. 
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Figure 2: Constraint from SNela observation in f2 rc -f2 m plane. Contours are for 95 % 
(dotted line) and 99 % (dashed line) C.L. constraints respectively. In this figure, we assume 
a flat universe. 
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